Cooperative Packet Routing using Mutual Information 

Accumulation 

Yanpei Liu, Student Member, IEEE, Jing Yang, Member, IEEE, Stark C. Draper, Member, IEEE 



Abstract — We consider the resource allocation problem in 
cooperative wireless networks wherein nodes perform mutual 
information accumulation. We consider a unicast setting and 
arbitrary arrival processes at the source node. Source arrivals 
can be broken down into numerous packets to better exploit 
the spatial and temporal diversity of the routes available in the 
network. We devise a linear-program-based algorithm which al- 
locates network resource to meet a certain transmission objective. 
Given a network, a source with multiple arriving packets and a 
destination, our algorithm generates a policy that regulates which 
nodes should participate in transmitting which packets, when 
and with what resource. By routing different packets through 
different nodes the policy exploits spatial route diversity, and by 
sequencing packet transmissions along the same route it exploits 
temporal route diversity. 

Index Terms — Wireless relay networks, cooperative communi- 
cation, mutual information accumulation, rateless codes. 



I. Introduction 

In multi-hop cooperative relay networks, relay nodes co- 
operate with each other to deliver data from network source 
to network sink. The benefits of cooperation include im- 
provements in energy efficiency, in robustness to fading and 
interference, and reductions in the likelihood of the loss 
of connectivity [1], [2]. Cooperation arises naturally in the 
wireless medium due to its broadcast nature. When two nodes 
are communicating, neighboring nodes can listen in on the 
transmission and thus are in a position to help. As an example, 
nodes further from the source and closer to the destination can 
overhear earlier nodes' transmissions. This gives them a head 
start on decoding the message and speeds up message delivery 
in contrast to conventional multi-hop transmission. 

In order to be able to exploit the broadcast nature of 
wireless transmission, we assume nodes are equipped with the 
following abilities: 

• Nodes can listen in on all transmissions. 

• Nodes store all observations and slowly accumulate noisy 
observations until they can decode the message. 

To simplify our model we assume transmitters operate on 
orthogonal channels and thus there is no interference. The 
first point simply means nodes continually monitor the full 
bandwidth of the system. The second point means that nodes 
have large memories and thus can store all observations 
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until they become useful. This is what we mean by mutual 
information accumulation. Mutual information accumulation 
can be realized through the use of rateless (or "fountain") 
codes [1], [3]-[7]. 

Enabling nodes with this physical layer (PHY) ability 
changes the routing and resource allocation problem from that 
traditionally considered (e.g., in backpressure [8]-[10] or in 
network coding [11]). Simpler PHY layer modications such 
as energy accumulation using space-time or repetition coding 
have been considered in [12]-[14]. The difference between 
energy accumulation and mutual information accumulation 
can be easily understood from the following example [1]. 
Consider binary signaling over a pair of independent erasure 
channels each having erasure probability pe from two relays 
to a single receiver. If the two relays use repetition coding, 
corresponding to energy accumulation, then each symbol will 
be erased with probability pi . Therefore, 1 — symbols are 
successfully received on average per transmission. Instead, if 
we use different codes, the transmissions are independent and 
on average 2(1 — Pe) novel parity symbols are received per 
transmission. 

In contrast to previous works that consider the routing of 
a single packet [1], [15], in this work we consider general 
arrival processes and the routing of multiple packets. This 
generalization is important to consider. With multiple packets 
on the go there are extra spatial and temporal degrees of 
freedom to consider that do not arise in the single packet 
case. As we show in this paper, by using our algorithm 
designed for multiple packets transmission, we can exploit the 
route diversity of the network in both spatial and temporal 
aspects. Different packets may be routed along partially (or 
wholly) disjoint routes, thereby exploiting spatial or route 
diversity. This is a desired feature when packets compete for 
resource; it may be desirable to route a group of packets 
through one link and the rest through another to balance the 
traffic and resource consumption. Depending on the arrival 
process different packets may be sent in sequence along the 
same route, thereby exploiting temporal diversity. This usually 
happens when routes choices are limited or packets arrive at 
different times. Essentially, our algorithm exploits available 
resource within the network and distributes it efficiently. 

The main contribution of our work is the design of a joint 
routing and resource allocation algorithm for multiple packets 
transmission from the source to the destination in a given net- 
work. Our work is a natural extension of the work in [1] where 
only single packet transmission is considered. The contribution 
of our work is twofold. We first present a formulation of 
the routing and resource allocation problem. The formulation 
considers various forms of energy and bandwidth constraints 
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and is expressed in the form of linear program (LP). Second, 
we design a centralized algorithm which solves the routing and 
resource allocation jointly by solving a sequence of LPs. Each 
LP solves for the optimal resource allocation given a route 
decision of all packets. The resource allocation result is then 
used to update the route decision and the method proceeds 
iteratively. 

Our work differs from previous cooperative routing and 
resource allocation works in the sense that we consider both 
PHY layer and routing (rather than just the PHY layer in 
simple networks or both in simplified settings). It differs 
from the current state-of-the-art cooperative routing algorithms 
(e.g., opportunistic multi-hop routing [2], backpressure-based 
algorithm [8]-[10], and network-coding-based algorithm [16], 
[17]) in that it involves a modified PHY layer with an extra 
feature - mutual information accumulation which can be easily 
implemented using rateless codes. Meanwhile, there are many 
existing works in cooperative routing using this feature. In [3], 
mutual information is considered for single relay networks. 
In [18], mutual information is also considered, but without 
a consideration on resource allocation. Routing and resource 
allocation with mutual information accumulation is considered 
in [1], [15], but they only consider single packet routing. 
To the best of our knowledge, there has been little prior 
work investigating routing and resource allocation with mutual 
information accumulation for multiple packets. 

The rest of the paper is organized as follows. We present 
the system model and problem formulation in Sec. |ll| The 
centralized algorithm is developed and discussed in Sec. [111 
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We provide detailed numerical results in Sec. IV Discussion 
and conclusion are in Sec. (V] 

II. System Model 

We consider a system with L nodes: the source, always 
labeled 1, the destination, always labeled L and L — 2 relay 
nodes. Suppose that K large data files (not necessarily of the 
same size) arrive at the source at a set of specified times 
Ti, r2, . . . , r^. Upon arrival each file is subdivided into a 
number of packets. In all there are > i^T packets. The 
indices of the packets are ordered according to the arrival 
times of the corresponding data files. We consider each of 
these packets as a separate commodity from a routing point 
of view. The network's goal is to deliver N such packets from 
the source to the destination under various channel constraints. 
Each relay may participate in transmitting a subset of these 
packets or remain silent. To simplify the analysis we assume 
that the only significant power expenditure for each node lies 
in transmission. 

A. Problem parameterization 

We first introduce the definition of a "decoding event". 

Definition 1: A decoding event occurs when either a node 
decodes a packet, or a packet is made available at the source. 

We assume a decode-and-forward relaying strategy. We 
denote as the time at which node i decodes packet c, 
c G {1,2,...,A/'}. In particular, each Tf corresponds to the 
arrival time of the corresponding data file I which packet c is 
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Fig. 1. A sample decoding order for = 3, L = 3 and M = NL. 



subdivided from, i.e., Tf = r^. We next present the definition 
of a "decoding order" - a key component in our algorithm. 

Definition 2: A decoding order T, containing Tf and 
T£ for all c e {1, 2, . . . , TV}, is a subset of the set 
{T^^, . . . , Tf^ . . . , T^} sorted in increasing order. 

It should be noted that a decoding order T only specifies the 
order of its decoding events. It does not specify when those 
events occur (or the value of s). 

We let T denote the decoding order and let M = | T\ denote 
the number of decoding events inside this decoding order T. 
Since a node cannot transmit a packet until it has decoded it, 
the order of decoding events puts constraints on the resource 
allocated to nodes for each packet. We further note that a 
decoding order can contain up to NL decoding events, i.e., 
M < NL. Then the label pair (z, c) of uniquely determines 
its position 5, s G {1,2,..., M}. In other words, there exists 
a unique mapping / which maps (z,c) to s and we further 
denote this mapping as T[ = Tj(^ c) = ^s- 

With a decoding order T, we define the "inter-decoding 
intervals" for its decoding events. 

Definition 3: Given a T, the inter-decoding interval A5 is 
the time interval between Tg, Tg-i G T and is defined as 

As = Tg — Ts_i, 

where we denote time as Tq. 

The entire message transmission can be thought of as con- 
sisting M intervals. The 5th interval is of duration A 5 and is 
characterized by the fact that at the end of this A5, a decoding 
event associated with T5 must occur. A typical decoding order 
is shown in Fig. [T] 

We assume the ith node operates at a fixed power spectral 
density (PSD) denoted as Pi (joules/sec/Hz), uniform across 
its transmission band. We assume the channel between any 
pair of nodes is block- fading and frequency non- selective. The 
channel gain between node i and node j is denoted as hij. 
Under these assumptions, the spectral efficiency (bits/sec/Hz) 
between node i and node j can be expressed as [1], [19] 
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where Nq/2 denotes the PSD of the white noise process. We 
further denote the time-bandwidth product allocated to node i 
to transmit packet c during a given time interval Ag as 
(sec-Hz). Then the information flow from node i to node j 
during this interval is ^Cij bits. 
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B. Problem constraints 

For a given decoding order we find the resource allocation 
minimizing an objective function subject to the following 
constraints: 

1) A5 > for all s. 

2) > for all s, i and c. 

3) Arrival process constraint; the source cannot transmit 
packet c until its corresponding data file has reached the 
source. 

4) Decoding constraint; node j must decode packet c at the 
associated timing point T5. 

5) Constraint(s) on energy and bandwidth. 

Arrival process constraint puts a constraint on value of Tf 's 
(thus the time intervals between them). Each Tf corresponds 
to the arrival time of the corresponding data file £ which packet 
c is subdivided from, i.e., Tf = r^. We describe constraints 
4) and 5) in more details below. 

We assume that nodes use codes that are ideal in the sense 
that they fully capture this potential flow, working at the 
Shannon limit at any rate. We also assume that distinct packets 
can be simultaneously transmitted from any node with no 
interference between packets and receivers have perfect knowl- 
edge to distinguish them. Nodes are further designed to use 
independently generated codes. This design leads to another 
assumption that a receiver can combine information flows from 
two or more transmitters for each packet without any rate loss. 
As noted in [1], the use of independently-generated codes 
is crucial for the mutual information accumulation process. 
If the transmitters used the same code, the receiver would 
get multiple looks at each codeword symbol. This is "energy 
accumulation." By looking at different codes the receiver 
accumulates mutual information rather than energy. 

With mutual information accumulation, the decoding con- 
straint imposes a constraint on some timing points. If Tj^, 
j 7^ 1 is in the decoding order, node j must be able to decode 
packet c at this time. This constraint is formally expressed as 

E E ^lsC^^3 > for an c, j / 1, (2) 

i:/(i,c)</(j,c)s=/(i,c) + l 

where 5^ is the size of packet c. Recall that Ci^j is the spectral 
efficiency (bits/sec/Hz) of the channel connecting node i to 
node j. Equation ^ says that in order for node j to decode 
packet c, the total accumulated information at node j for 
packet c must exceed bits by the time TJ. The non-ideal 
nature of existing implementations of rateless codes can be 
handled by incorporating an overhead factor (1 + e) into the 
right-hand side of ([2]). Also, the relation ^ expresses the 
constraint that node i can transmit packet c to node j only 
after it has decoded packet c. 

Constraints on energy and bandwidth can either be system- 
wide constraints or be imposed on a node by node basis. We 
state various possible constraints in the followings. 

1) Per-node bandwidth constraint. If node i is assigned 
with bandwidth Wi, its resource allocated during given 



A 5 must satisfy the following 

N 

Y,Als<^sW,, for an 5. (3) 

c=l 

2) Sum bandwidth constraint: If the total bandwidth Wt 
is allocated, a sum bandwidth constraint applies across 
all nodes during any given A5. This can be formally 
expressed as 

L N 

E E ^ ^sWt, for an s. (4) 

i=l c=l 

3) Per-node energy constraint. If node i with transmission 
power Pi is assigned with energy budget Ei, we can 
express the per-node energy constraint as 

N M 

E E ^isP^ ^ 

c=l s=/(i,c) + l 

4) Sum energy constraint. If a sum energy Et is assigned 
for all nodes, a sum energy constraint can be applied as 

L N M 

E ^lsPi<ET. (6) 

i=l c=l s=/(i,c)+l 

In the next section we develop a centralized algorithm based 
on the LP framework. 

III. Centralized Algorithm 

The LP framework can take many objective functions. If 
there is only one data file (thus T^^Tf , . . . = ri), one 
natural objective function to use is the total transmission time 

M 

Alternatively, if K data files arrive at different times (thus 
packets are not available to the source at the same time), 
instead of minimizing total transmission time, one might want 
to minimize the average transmission time given by 

Ta = ^f2Te, (8) 

£=1 

where is the duration data file i is being transmitted. This 
objective function is appropriate in the sense that it does not 
penalize the arrival time. 

Other linear programming frameworks are also possible. 
For example, one may wish to minimize the total energy 
expenditure given by 

L N M 

EE E (9) 

i=l c=l s=/(i,c) + l 

subject to total transmission time or average transmission time 
constraint. 
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A. Characteristics of the problem 

We now study the properties of routing and resource alloca- 
tion under different constraints with total transmission time ^ 
as our objective function. First, under the per-node bandwidth 
constraint, we have the following lemma. 

Lemma 1: Under a given decoding order, suppose and 
^ are the solution to the problem under the per-node 
bandwidth constraint ([Sj. Then for each 5, there exists a node 
i such that the inequality in ([3| becomes equality. 

Proof: We prove by contradiction. Suppose for some s 
inequality is strict for all i, i.e., 

N 



^Al^ </\sWi, foralH. 



Then we can scale down A5, yielding a smaller objective 
value. This contradicts the assumption that A5 is the optimal 
solution to the LP. ■ 
This proposition suggests that the A5 can be calculated as 



max 



Z^c=l ^i,s 



Under the sum bandwidth constraint, we first have the 
following lemma. 

Lemma 2: Under a given decoding order, suppose A5 and 
Al g are the solution to the problem under the sum bandwidth 
constraint (|4]). Then the inequality in ^ must be equality for 
all s. 

Proof: We prove this by contradiction. Suppose for some 
s the equality does not hold, i.e.. 
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AsWs < AgVFr, for some s. 



Then for these s we can scale down the corresponding A 5 by 
^ while increase Wg to Wt- We therefore obtain a solution 
which has smaller objective value. ■ 

This lemma leads to the following important theorem. 

Theorem 3: Under the sum bandwidth constraint and a 
given decoding order, if Pi = P for all i then the solution 
that minimizes the objective in ^ also minimizes the sum 
energy. 

Proof: The total energy expenditure of the entire system 

is 

L N M 

EE E ^isP 

i=l c=l s=/(i,c) + l 
L N M 

i—l c—1 s=l 
M 

= ^AsWtP 
= TtWrP^ 

The first equality holds because a node cannot transmit a given 
packet (resource allocated to it is zero) before it decodes the 
packet. The second equality follows from Lemma [2] Since the 



objective Tt is proportional to the energy used, minimizing 
one minimizes the other. ■ 
This theorem tells us that in this setting there is no trade 
off between total transmission time and energy. The minimum 
transmission time policy is identical to the minimum energy 
policy. 

Theorem 4: Under the sum bandwidth constraint and a 
given decoding order, if the minimum transmission time for 
transmitting one packet is T, then the minimum transmission 
time for routing N packets of same size is NT. 

Proof: Let 71 denote the optimal set of routes on which 
transmitting a packet takes T amount of time. These routes 
are equivalent in the sense that transmitting a packet on any 
of them takes T amount of time, consuming TWt a amount of 
total resource. Then transmitting TV packets of the same size 
on 7Z can be finished within NT amount of time by allocating 
each packet amount of total resource. 

We now show that this is the best the system can do. 
Suppose otherwise, i.e., not all packets are routed through 
TZ and the total transmission time is smaller than NT. Then 
there exists at least one packet traveling on a different set 
of equivalent routes 1Z' ^ IZ. Since the transmission time 
is smaller than NT, each packets traveling on IZ should be 
given more than amount of resource. This implies that 

the packet on IZ' shares less than amount of resource 

whereas achieves a total transmission time less than NT. 
This contradicts our assumption that IZ is the optimal set of 
equivalent routes. ■ 

This theorem suggests that in sum bandwidth constraint, the 
performance may not be improved by dividing a big data file 
into smaller packets. We will illustrate this in simulation. 

B. Optimizing the Decoding Order 

In this section we show an important theorem that moti- 
vates the centralized algorithm. The theorem tells us how to 
manipulate the decoding order based on the solution of the LP 
problem. Denote as 7^ a LP framework with a chosen objective 
and constraints being considered. Given any decoding order of 
length M, define 



Ai,...,AM,A{,i,...,Af, 



M 



to be the solution to the LP with optimal objective value Topt 
on this decoding order. We then have the following theorem. 

Theorem 5: If A^ = for some m and we swap the 
positions of and T^_i in the decoding order, then the 
objective value T*^^ obtained with this swapped decoding 
order satisfies T*^^ < Topt- 

Proof: We prove the theorem by showing that when 
swapped decoding order is used, the original solution x with 
optimal objective value Topt is still feasible under new decod- 
ing order. To show x is feasible under the new decoding order, 
we show the decoding constraint of the swapped decoding 
order is satisfied by this solution. 

If a decoding event associated with T^-i has its decoding 
constraint satisfied at this time, its decoding constraint is 
certainly satisfied at a later time Tm- If a decoding event 
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associated with has its decoding constraint satisfied at this 
time, it actually has the constraint satisfied at an earlier time 
Tm-i since = for all i and c. Therefore, solution x 

with optimal objective value Topt is feasible under swapped 
decoding order ■ 
The idea behind Theorem |5] is based on the following 
observation. A solution to the LP with = indicates 
that the event associated with Tm takes place at exactly the 
same time with the previous event T^-i, or actually occurs 
before it. Therefore, swapping the position of and Tm-i 
typically gives a decrease in the objective value once the LP 
is solved under swapped decoding order. In the case that T[ is 
swapped with T£ for i not equal to L, the decoding event T[ 
is excluded from the new decoding order by our algorithm. 

C. Centralized algorithm 

We design and implement a centralized resource allocation 
and routing algorithm iterating between two sub-problems: 

1) For a given decoding order, a resource allocation scheme 
is determined by solving a given LP problem V. Solving 
the optimization problem with these constraints gives 
time allocation (thus Tg) for all s and the resource 
allocation ^ for slU i, c and s. 

2) With a given resource allocation scheme, we update the 
decoding order. 

In the following we discuss the algorithm in more details. 

1) Initialization: We first consider the decoding order ini- 
tialization as if there were only one packet to be routed. 
We design a heuristic algorithm similar to the greedy filling 
algorithm in [12] based on network state. For single packet 
transmission, the decoding order is essentially the order of 
nodes which successfully decode that packet and can come 
online to broadcast it. Therefore, with a set of already decoded 
nodes, the next node to decode is most likely the one which is 
able to benefit the most from the already decoded nodes. We 
thus choose the next node to be the one that has maximum sum 
spectral efficiency between itself and all previously decoded 
nodes. The algorithm pseudocode is given in Algorithm [T] 

Algorithm 1 Initialize the decoding order 

1: T^{Tl} 

2: while length{T) < L — 1 do 
3: j ^ arg maxfe^r,fe/L EieT 
4: T^{T,Tl} 
5: end while 
6: T^{T,n} 



In multiple packets case, since we have no initial knowledge 
about which nodes should decode which packets and in which 
order, we believe it is difficult to design an effective decoding 
order initialization algorithm. We therefore initialize the de- 
coding order based purely on the one with single transmitted 
packet. Suppose in single packet case the following decoding 
order is found 

r=lT^,Ti,---,Tl], 



the decoding order for multi-packet case is simply chosen to 
be 

r=|T^...,T,^,7;\...,7;^,...,Ti^|. 

2) Decoding order updates: Based on the solution to the 
LP, the algorithm first searches for all those m's such that 
A^ = 0. Then for those m, it checks for the situation where 
T[ is followed by T£ for some i ^ L.lf it finds one, it drops 
the corresponding T[ and reruns the LP. If it does not find one, 
it swaps the order of Tm-i and Tm- If the swapped decoding 
events are all the same as the ones in the previous iteration, 
and if there was no decoding event dropped in the previous 
iteration, the algorithm terminates. We present the sketch of 
the algorithm in Algorithm |2] 

Algorithm 2 The routing and resource allocation algorithm 

1: Choose a LP 7^ 

2: Initialize a decoding order T 

3: Do Again ^ true 

4: while DoAgain do 

5: result ^ LP(T, V) 

6: index ^ SearchForZeroDeltas(result) 

7: ToDrop ^ SearchForNodeDecodeAfterDestina- 
tion(index) 

8: if ToDrop 7^ then 

9: DoAgain ^ true 
10: drop the decoding event 

11: else if drop = and index = index in last iteration and 

no event is dropped in last iteration then 
12: DoAgain ^ false 
13: break 
14: else 

15: DoAgain ^ true 
16: SwapEvents(result) 
17: end if 
18: end while 



Because of the exponential number of orderings we expect 
the problem of finding the optimal decoding order to be NP- 
hard. The sub-optimality of our heuristic algorithm comes 
from the fact that the excluded decoding events may actually 
be helpful. Without a mechanism to "re-introduce" the ex- 
cluded events our algorithm is not expected to achieve global 
minimum for all networks. 

3) Characteristics of final route: Our algorithm has the 
property that different packets may take different paths toward 
the destination under the per-node resource constraint. This 
exploits the spatial route diversity. For example, the algorithm 
may schedule the first 3 packets to pass through a set of relay 
nodes and schedule the next 3 packets to pass through yet 
another set of relay nodes. The intuition is that instead of 
choking an existing route and overloading the already busy 
nodes, the algorithm balances the traffic by invoking other 
free nodes to fully utilize the network. 

4) Challenges: The challenges we face in designing this 
multi-packets routing and resource allocation algorithm are 
more difficult than the ones the authors encountered in [1]. 
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Fig. 2. A diamond network example. Note that transmitting a single packet 
on route 1 — 3 — 4 or 1 — 2 — 4 takes the same amount of time. 



Essentially, we are trying to sort out a correct order for 
both nodes and packets. Our algorithm may be stuck in local 
solutions without the ability to escape. However, as we show in 
simulation, for smaller networks, although it may not always 
find the optimal solution, the local solution it achieves is still 
promising. 

IV. Numerical Results 

In this section we present detailed simulation results for 
the algorithm. These results exemplify the basic properties of 
using mutual information accumulation in cooperative com- 
munication. 



A. A simple network 

To better illustrate how our algorithm exploits the spatial 
route diversity, we first consider a simple four-node example 
in which a data file of size 20 is equally packetized into two 
packets of size 10. Consider the diamond network in Fig. |2] 
Links with nonzero spectral efficiency (bits/sec/Hz) are shown. 
Each node is allocated with 1 unit of bandwidth and the source 
(labeled as 1) transmits these two packets to the destination 
(labeled as 4). 

Running our algorithm returns the following routing and 
resource allocation scheme. Starting from time 0, node 1 
broadcasts packet 1 for 5 units of time. By the end of time 
5, node 2 decodes packet 1 and node 3 accumulates half of 
packet 1. Starting from time 5, node 2 broadcasts packet 1 
to node 4 while node 1 broadcasts packet 2. By the end of 
time 15, node 4 decodes packet 1 and node 3 decodes packet 
2. Note that by the end of time 10, node 2 decodes packet 
2. Then starting from time 15, node 2 and 3 transmit packet 

2 to node 4 at the same time and node 4 decodes packet 2 
after 3.33 units of time. The total transmission time is thus 
5 + 10 + 3.33 = 18. 33 units of time. In this example, node 

3 never decode packet 1 thus packet 1 does not pass through 
route 1 — 3 — 4. Whereas node 3 decodes packet 2 and it 
cooperates with node 2 to deliver packet 2 to destination. 

It is interesting to see that if two packets were routed 
through the same path, the total transmission time would be 
2 X 15 = 30. If the data file itself is routed without packetizing, 
the total transmission time would be 30 no matter which 
route (1 — 3 — 4 or 1 — 2 — 4) it takes. This suggests that 
possible improvement can be achieved if we divide data files 
into smaller packets. We further illustrate this property in the 
following section. 



Empirical CDF 

■"l ' ■ : 

0.9- y^Y 

// I / 

0.8- 77 I / 

0.7- # /; / \ / 

0.6- // / ; / ; / 

; // ' ; / : / 



0.3- 
0.2- 



N 


= 1 pacl<et 


N 


= 2 pacl<ets 


N 


= 3 pacl<ets 


---N 


= 3 Dijkstra with MIA 


N 


= 3 Dijkstra without MIA 



°0 5 10 15 20 25 30 

Total Transmission Time 

Fig. 3. Cumulative distribution of total transmission time under N = 1, 
2, and 3. While networks randomly generated within the unit square have 
constant number of uniformly distributed nodes, the Dijkstra routes on these 
networks may contain different number of nodes. This contributes to the non- 
smoothness in the Dijkstra algorithm simulation plot. 

B. Simulation on general networks 

We now study the performance of our algorithm on more 
general networks. For better illustration, we choose bandwidth 
as our single resource element being allocated. 

We simulate the performance of our algorithm on 10 
randomly scattered wireless nodes within a unit square. The 
source, denoted as 1, is always located at [0.2,0.2] and the 
destination, denoted as L, is always located at [0.8,0.8]. The 
remaining 8 nodes are randomly placed according to the 
uniform distribution in the unit square. The channel gain hi^j is 
deterministically related to the Euclidean distance di^j between 
node i and j as hi^j = d~^. The spectral efficiency between 
node i and j is given by 

1) Per-node bandwidth constraint: We first consider the 
system with the per-node bandwidth constraint. We consider 
a system with single arrival of a data file of size 20 which 
is equally packetized into N packets. Each packet is of size 
20/7V. We also let Pi = P = 1, Wi = 1 and No/2 = 1. 
We simulate networks to estimate the cumulative distribution 
function (CDF) of the total transmission time when we either 
do not subdivide the data file (i.e., the same problem as con- 
sidered in [1]) or divide it into 2 or 3 packets for transmission. 
To further illustrate the benefits of using mutual information 
accumulation, we also plot the transmission performance on 
the route obtained using Dijkstra's algorithm [20] with and 
without mutual information accumulation in Fig. |3] 

We see that dividing a data file into packets improves the 
performance by reducing the total transmission time. If we 
keep dividing it into more and more packets, we expect that the 
performance eventually converge as if the data file is treated 
as fluid. Comparing with Dijkstra's algorithm, we also notice 
that the decrease in total transmission time comes from two 
causes; the use of mutual information accumulation and the 
ability to exploit other available routes and resource within a 
network, not merely confined in a pre-determined Dijkstra's 
route. This illustrates how our algorithm exploits the spatial 
route diversity. 
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Fig. 4. A randomly generated network. Transmitting a single data file of 
size 20 on this network takes 7.1576 units of time. 
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Fig. 5. Average transmission time under different inter-arrival times. 

2 ) Per-node bandwidth constraint with arrival process: We 
now consider a deterministic arrival process under the per- 
node bandwidth constraint. Consider a randomly generated 
network shown in Fig. |4] We now route K = 3 data files 
with inter- arrival time between two consecutive files, each is 
of size 20, from the source to the destination. We also let 
Pi = P = I, Wi = I and Nq/2 1. Instead of using 
total transmission time as our objective, we try to minimize 
the average transmission time given in ([5]). We simulate the 
average transmission time required when inter-arrival time 
(7"2 — Ti or Ta — T2) varies from to 15. 

The simulation shows the exploitation of the temporal route 
diversity. For this given network, transmitting a single data file 
of size 20 requires 7.1576 units of time. If the inter-arrival 
time is larger than this number, then transmitting 3 arriving 
data files means repeating the single file routing and resource 
allocation policy three times. Thus the average transmission 
time is the same as transmission time for a single data file. 
When the inter-arrival time is small, files compete with each 
other for resource, thus the average transmission time increases 
as shown in Fig. [5] 

We now divide each data file equally into 2 packets of size 
10 (thus N = 2K = 6). Since under the per-node resource 



Empirical CDF 




Total Transmission Time 



Fig. 6. Simulation of the sum bandwidth constraint and comparison with 
the per-node bandwidth constraint. 

constraint dividing a big data file into packets better utilizes 
network resource, we have a performance gain shown in Fig.|5] 
Note that the inter-arrival time is between every two packets 
of size 10. 

3) Sum bandwidth constraint: We next simulate the per- 
formance of the system with the sum bandwidth constraint. 
We consider a single data file of size 20 which is equally 
divided into N packets. We let Wr = 10, = P = 1, 
No/ 2 = 1. Note that the total bandwidth is the sum of 
individual bandwidth Wi in the per-node bandwidth constraint 
simulation setting. 

The simulation result is shown in Fig. [6] Interestingly, we 
note that under the sum bandwidth constraint, dividing a data 
file into more and more packets would not help to decrease 
the total transmission time. This is the result of Theorem |4} 
since each packet is also ^ of the size of the data file, the 
total transmission time would not change. 

We also notice that in Fig. [6j the performance under the 
sum bandwidth constraint is better than that under the per- 
node bandwidth constraint. The reason is that in allocating 
each node a fixed amount of bandwidth as in the per-node 
constraint, the nodes which cannot transmit waste the resource 
assigned to them. In the sum bandwidth constraint, the network 
gathers all available resource and distributes it to serve the 
most needed nodes; the nodes which cannot transmit are not 
assigned with resource. 

V. Summary and Future Work 

In this paper we study the problem of cooperative commu- 
nication using mutual information accumulation. We design 
and prototype a routing and resource allocation algorithm 
based on solving a LP based problem iterating between two 
subproblems; finding the best decoding order and finding the 
best resource allocation given a decoding order. Our algorithm 
exploits both the spatial and temporal route diversity and 
suggests that performance can be improved by dividing a 
big data file into several small packets. Our algorithm also 
supports packets with different arrival times. 

In some scenarios, it is not possible to centralize the routing 
and resource allocation scheme. Limitations on centralized 
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algorithm are constraining in large scale networks where ag- 
gregating the channel state information can incur unacceptable 
overheads. The centralized algorithm is also not desired in 
temporally varying networks where central controller often has 
trouble in updating the estimates of channel state information 
promptly. Developing a distributed algorithm would be an 
important direction for future research. 

Conventional wireless nodes have the ability to adjust the 
transmitting power. With individual transmitting power as a 
varying variable, it is therefore natural to consider power 
allocation optimization. We note that individual transmitting 
power affect the spectral efficiency according to ([T]). Thus 
maintaining both the time-bandwidth product and the transmit- 
ting power as optimization variables makes the decoding con- 
straint nonlinear. Understanding how power allocation would 
affect the transmission performance is an interesting topic for 
investigation. 

Consider a cooperative broadcast problem in which the 
objective of the network is to deliver packets from the source 
to all other nodes. Since every node is a destination, we do not 
need to drop a decoding event. This preserves some optimality 
of our algorithm. However, it is still unknown that if our 
algorithm is the most effective one. It would be interesting 
to investigate possible extensions of our algorithm tailored for 
this problem. 
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